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Abstract. In general relativity all forms of energy contribute to gravity and not only just ordinary matter as in 
Newtonian Physics. This fact can be seen in the modified hydrostatic equilibrium equation for relativistic stars 
pervaded by magnetic (B) fields. It has an additional term coupled to the matter part as well as an anisotropic 
term which is purely of magnetic origin. That additional term coming from the pressure changed by the radial 
component of the diagonal electromagnetic field tensor, weakens the gravitational force when B is strong enough 
and can even produce an unexpected change in the attractive nature of the force by reversing its sign. In an extreme 
case, this new general relativistic (GR) effect can even trigger an instability in the star as a consequence of the 
sudden reversal of the hydrostatic pressure gradient. We suggest here that this GR effect may be the possible 
central engine driving the transient giant outbursts observed in Soft Gamma- ray Repeaters (SGRs). In small 
regions of the neutron star (NS), strong magnetic condensation can take place. Beyond a critical limit, these 
highly magnetised bubbles may explode releasing the trapped energy as a burst of 7-rays of ~ 10 36-40 erg. 
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1. Introduction 

A class of neutron stars seem to have strong surface 
magnetic fields of the order of 10 15 Gauss. They are re- 
ferred to as the Magnetars. These highly magnetised stars 
have been speculated to be the progenitors of some pe- 
culiar gamma ray bursts events, the so-called repeaters. 
Two kind of (candidates) compact stars are believed to 
have very strong surface magnetic field: the Soft Gamma- 
ray Repeaters (SGRs) and the Anomalous X-ray Pulsars 
(AXPs: Mereghetti & Stella 1995). SGRs are objects that 
repeatedly emit bursts of gamma rays with energies of the 
order of io 36 ~ 40 erg, in addition to persistent X-rays. The 
SGRs are also a small enigmatic class of hard X-ray tran- 
sient sources (Marsden et al. 2001). AXPs show persistent 
X-ray emission, modulated at a stable, slowly lengthening 
period. Contrary to the standard binary X-ray pulsars, 
they show no evidence for a companion star. 
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Stars with strong surface magnetic fields have still 
stronger field inside. Influence of strong magnetic fields 
on degenerate nuclear matter has been studied by Lai & 
Shapiro (1991). They employed the scalar virial theorem 
2T + W + 3U + M =0, withT, W, II and Al, being the ro- 
tational kinetic energy, the gravitational potential energy, 
the internal energy and the magnetic energy respectively. 
With this approach, they found out that a surface field of 
10 13 Gauss can readily imply that an inner magnetic field 
to be as high as 10 18 Gauss. Alternately, this field limit 
can easily be estimated if we consider the magnetic energy 
approaching the gravitational energy to produce a signifi- 
cant deviation from the normal balance between pressure 
and gravity. This implies an extreme upper limit for the 
average field B, with (S 2 /87r)(47ri? 3 /3) - GM 2 /R so that 
B Um ~ 1.31O 8 (M/M )/(i?/i? Q ) 2 (Mestel, 1999). For a 
neutron star with R/Rq ~ 10~ 5 , Bu m ~ 10 18 Gauss. In 
this work we discuss the effect on the gravitational force, 
of strong magnetic fields of the order of the previously 
quoted limit. 
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Unlike its Newtonian counterpart, in the GR hydro- 
static equilibrium equation there appears a new term that 
depends not only on the pressure of the gas but also on 
the electromagnetic field strength in the radial direction. 
When the magnetic field is very strong and approaches the 
critical limit, it "softens" the gravitational force. In the ex- 
treme case, when the threshold strength (10 18 Gauss) is 
reached, the sign of the radial pressure gradient re- 
verses, thus indicating that the overall effect has become 
repulsive. An instability must then appear, and a crite- 
rion for it can readily be set. The appearance of a positive 
pressure gradient in the radial direction, i.e., > 0, 
means that the NS is unstable. This may be interpreted 
as if in the interior of the magnetar the matter struggles 
to expand as the field reaches (locally) such large values. 

The electromagnetic field can be treated as a fluid and 
there exists room for the occurrence of magnetic field con- 
densation in a local volume in the NS interior. The B limit 
is called for here to study highly magnetised matter bear- 
ing in mind a stable stellar configuration. So, it is not sur- 
prising that in order to reach the instability one needs to 
be nearly at this critical field limit. The radial instability 
discussed here, and expected to occur in the interior of a 
magnetar, is actually produced by the spatial anisotropy of 
the magnetic stress-energy tensor. The anisotropic stress 
tensor (i.e., the fact that the radial component, p r , differs 
from the angular components, pe = P4, = Pt or, the trans- 
verse component), introduces a repulsive term, purely of 
magnetic origin, in the hydrostatic equilibrium equation. 
Besides this term, the stellar gravitational pull is reduced 
because it depends directly on the radial component of the 
stress-energy momentum tensor which becomes smaller in 
the presence of a very strong magnetic field. Both effects 
conspire to make positive the "effective" gradient of the 
radial pressure as to induce the instability. 

A similar effect has been discussed recently by Ray et 
al. (2003, 2004) with regard to the appearance of a net 
electric charge in a stellar structure. Therein, it is shown 
that there exists a critical limit to the electric field, beyond 
which the star is unstable. Here we make use of this novel 
phenomenon, the counterbalance of the gravitational pull 
excerpted by the compact object in critically magnetised 
interiors of neutron stars. The attentive reader must re- 
alize that not everywhere in the interior of such stars the 
field is so high as ~ 10 18 Gauss, except only in small re- 
gions (the bubble) where it can be raised by the effect 
of a strong magnetisation, as discussed for example by 
Perez Martinez et al. (2003). This instability may trigger 
an abrupt release of energy similar to the giant flares of 
soft gamma-ray outbursts seen in SGRs. 

The mechanism for the emission of gamma rays from 
the SGRs is still unknown. Typical energy of the emitted 
flares in the SGRs is 10 36 erg/s, with special exception of 
a few sources emitting as high as 10 44 erg/s. The mag- 
netar model of Wheeler et al. (2000) makes use of the 
high magnetic field (10 15 to 10 17 Gauss) of a rapidly ro- 
tating neutron star as the possible source of cosmological 
Gamma-ray bursts (GRBs) releasing energy of ~ 10 54 erg. 



As a manifestation of the radial instability generated by 
the GR effect discussed here, we also suggest a model to 
the emission process of the SGRs, by the concentration 
of the magnetic field in certain volume of the star, which 
we may call it as a "magnetic bubble" . It is the fast ex- 
pansion (explosion) of this bubble, driven by the positive 
pressure gradient, which triggers the quiescent state to the 
super outbursts in SGRs, and as such, it may be a viable 
explanation of their central engine. 

Magnetic field evolution in compact stars have been 
studied in details by Thompson & Duncan (1993, 1995, 
1996 & 2001) and by Thompson, Lyutikov & Kulkarni 
(2002). Their main idea for the existence of high field in 
the neutron stars (magnetars) was, due to some dynamo 
mechanism (alpha - Omega), based on the fact that the 
neutron star convection is a transient phenomena and has 
extremely high Reynolds number. They showed that most 
of the magnetic energy of young pulsars reside in a convec- 
tive cell of diameter < 1 km. Lyutikov & Blandford (2002 
& 2003) modeled the cosmological GRBs as explosions 
of electromagnetic bubbles in a force free configuration, 
where the bubbles expand non-relativistically inside the 
star, and then become highly relativistic when they break 
from the surface. 

Here also we take such an idea as a working hypothesis 
of the formation of small volumes like bubbles, where the 
magnetic field is very high and nearly close to some critical 
field limit. The order of this critical magnetic field inside 
the magnetic bubble is found to be the same as that of the 
pressure and the mass of the inner sphere where the bubble 
is located. Also the not so frequent giant flares found in 
SGR 0526-66, with the famous March 5, 1979 event of 
energy release of 1.3 x 10 44 erg, and the SGR 1900+14, of 
the August 27, 1998 event of energy release of 6.8 x 10 43 
erg, can also be caused by this relativistic effect of the 
strong magnetic field on the gravitational force. We argue, 
conversely, that what the magnetic field drives in fact is 
the reversal of the gravitational force when a critical field 
strength is reached. 

We organise this article in the following way. In 
Section we present the formalism to describe the stellar 
matter under strong magnetic fields in the light of general 
relativistic formulation and present the modified Tolman- 
Oppenheimer-Volkoff (TOV) equation. In the Section |3J 
we explain the instability and estimate the critical value 
of B field as a function of the depth of the star. In Section^] 
we present the emission mechanism of the bursts in SGRs 
and give the estimates of the released energy and finally 
we draw our conclusions in Section 

2. General relativistic effects of magnetic fields on 
the structure equations 

In this section we study the effect of the magnetic field 
in the hydrostatic equilibrium equation. The magnetic 
field energy-momentum tensor introduces an anisotropy 
in the fluid pressure in the stellar interior. In our case we 
have an anisotropy in the sense that the radial component 
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of the pressure p r differs from the angular components 
Pe = P<f> = Pt (transverse pressure). We choose the mag- 
netic field produced by a magnetised surface to be a radial 
function which guarantees that the pressure components 
are still a function of the radial coordinate keeping the 
spherical symmetry. 

We proceed with the standard form of the metric 



ds 2 



e v dt 2 



e x dr 2 



r 2 {d9 2 



sin 



(1) 



The Einstein-Maxwell stress tensor will have con- 
tributions coming from the matter part and the magnetic 
fields, and will take the form: 

T£ = (p+e)u»u v ~pg»+±- f^-F» a F av + ^F a(3 F a ^j (2) 

where p is the pressure, e is the energy density (=pc 2 ) 
and u-s are the 4-velocity vectors. For the time compo- 
nent, one easily sees that u t = e u / 2 and hence u l u t — 1. 
Consequently, the other components (radial and spherical) 
of the four vector are absent. 

Maxwell's electromagnetic field equations can also be 
written as: 



Thus by using Eq.JSJ we have, 

r 2 B)' =0, 



(10) 



where the prime stands for the derivative with respect to 
the coordinate r. Such an equation can be integrated to 
give 

B(r)=q m /r 2 , r >0 , 

q m being an integration constant. 

The nonzero components of the stress-energy tensor 
T,f are 



T\ = (e + Pern), % = ~{p ~ Pern) 

and 



T<j = T* = -(p + Pem ) 



(11) 



where p em = B 2 /8ir is the energy density in the elec- 
tromagnetic field. The relevant components of Einstein- 
Maxwell field equations are 



f — n 



(3) 
(4) 



where F^- a and are the covariant derivative of the field 
strength tensor and the electric current 4-vector, respec- 
tively. Eq. (PJ gives the form of the vector potential (A^) 
as F^u — A v ^ — A^ u . So, the electric and magnetic fields 
in the rest frame are defined as 



Fa — Fpi/U , 



Bu 



1 

--e„ 



v pap 



(5) 
(G) 



C/ii/a/3 being the Levi-Civita antisymmetric unit tensor 
with £oi23 = \T~9- Now, considering the presence of a 
purely magnetic field, we can write 



F, v u v = 0. 



(X) 



This can also be thought of as a highly conducting 
medium of the star, where the term for the electrostatic 
charge will vanish. 

The electromagnetic field strength is chosen in the 
form 



Fa 



-F v g=B(r)r 2 sin (9, 



(8) 



where B(r) can be interpreted as the local magnetic field, 
the other components of F^ v being zero. 

Maxwell's equations imply just one equation, the t 
component 



W r F u = 4tt J* 



(9) 



,-A 
-X 



e " 1 
— (rA' - 1) + ^ = 8vr (e + p em ) 



e~" 1 
—2- [rv + 1) - — = 87r (p - p em ) 



(12) 
(13) 



The first of the Einstein's equations is used to deter- 
mine the metric coefficient e A , which is used to define the 
mass and charge inside a sphere of radius r. Namely, 



—A I _ 2m (r) <l 2 m {r) 

— 2 ' 

where 



(14) 



m(r) 



Anr 2 (e + p em ) dr + 



2r 



(15) 



The mass of the star is now due to the total contribu- 
tion of the energy density of the matter and the magnetic 
energy (B 2 /8n) density. Such a definition is suitable since 
it coincides with the mass of a localised object as observed 
from infinity. So, for an observer at infinity, the mass is 
given by 



R poo 

m - x = / Airr 2 (e + p em ) dr + / 4nr 2 (e + p em ) dr 

10 J R 

:R 4Trr 2 (e + p em )dr+ti^=m(R), (16) 



2R 



where R is the radius of the star. 

From Maxwell's equations we find q m (f) = B(r)r 2 — 
q m (0). If there are no singularities at r = the charge 
<7 m (0) is zero. Therefore, it follows the well known relations 
for radial magnetic field B(r) = q m /r 2 (as above). 
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Exact solutions for anisotropic stars were obtained 
by Dev & Gleiser (2002;2003), where they considered 
anisotropy in the pressure component for a spherically 
symmetric gravitationally bound object. Taking a radial 
and transverse component of the stress-energy tensor (p r 
and pt respectively), they arrived at the analogous TOV 
equation for anisotropic stars. 



dp r 
dr 



(p + p r ){M(r) +47rrV) 



2 (l- 



2M(r) 



+ ~(Pt ~Pr) 



(17) 



This equation explicitly shows the repulsive term com- 
ing from the difference between the transverse and ra- 
dial pressures (anisotropy effect), with the radial pressure 
changing the gravitational attraction. This is a purely rel- 
ativistic effect and it is of importance because of the ap- 
pearance of a radial magnetic field. To our knowledge, it 
has not been considered in previous studies of magnetars. 

Considering that T? = — pi for i = r,9 and <p an d 
using in the equilibrium equation, Eq. (|17|l . the stress- 
energy momentum tensor from Eq. (| 1 1 11 , and replacing the 
energy density p by the time component T\ = (e + p em ), 
the modified TOV for this case will be: 



dp r 
dr 



(P- 



M(r) + 47rr 3 (p - f£ 



£2_ 

27rr 



(18) 



where, the mass M{r) is the mass of the inner sphere, 
below the depth where the bubble is situated, and is inte- 
grated along the entire length, from the center of the star 
to the bubble. Here M(r) is the mass for an observer in 
the star and is related to the mass at infinity m(r) from 
Eq. (JT2J by M(r) = m(r) - q 2 m (r)j2r. 

This form for the modified hydrostatic equation equa- 
tion is essentially the same as has been recently done for 
the case of presence of charge inside a star (Ray et al. 
2003) if we substitute the radial pressure gradient by the 
total dp I ' dr using the expression for p r = (p — p em ) from 
Eq. JTTJ. 



dp r 

dr 
or. 



< 



(p + e) (M(r)+47rr 3 (p-f£)) 



D 2 



r 2 I i _ 2M(r) \ 2irr 



> 0. 



(19) 



This can be written in terms of the limiting magnetic 
field as 



\ 



2tt{p + €)- 



(M(r) + 4irr 3 p) 



r 1 - 



2M(r 



+ 7r(p + e)r 2 



(20) 



Magnetic field stronger than this critical field limit re- 
verses the sign of the radial pressure gradient, as shown 
in the Fig. (JTJ. 




4 6 8 

Radius [km] 



Fig. 1. The dp r /dr of the star as a function of the radius 
of the star. The region of the space where it faces a strong 
magnetic field, the dp r /dr becomes positive, and produces 
a burst. The maximum strength of the field is necessary 
where the depth of dp r /dr is the maximum. 



3. A radial hydrostatic instability created by 
strong magnetic fields 

Let's find out the order of magnitude of the magnetic field 
so as to produce an instability inside the star. Considering 
Eq. I|18fl . the negative sign of the radial pressure gradient 
needs to be maintained for the entire radius of the star 
for a stable configuration. The reversal of the force field 
occurs only for a large value of the magnetic field. Now, the 
first term on the right hand side of Eq. I|18fl . i-c, (p + e) is 
positive definite. So, the term which primarily controls the 
'change of sign' for dp r /dr is (M(r) + 4irr 3 (p — B 2 /8tt)). 
Also, the presence of the positive term B 2 /(2nr) helps in 
the process of making the dp r /dr positive. So the stability 
criterion is 



What is the order of the critical magnetic field required 
to equate both sides of the inequality (|20|l ? Taking pres- 
sure (p) to be in MeV/ fm 3 , relation (|20|l shows that the 
analog of B in that unit is (MeV ' j fm 3 ) 1 / 2 . Hence, 

1 (MeV/fm 3 ) 1 ' 2 ~ (1.6 x 10 33 ) 1/2 c..g.s. units 

~ 4 x 10 16 Gauss. (21) 

Working out relation H20|) we find that for the 
characteristic pressure inside a neutrons star p ~ 
10(MeV/ fm 3 ) ~ 10 35 dynes / cm 2 the critical field limit is 
- 10 2 (MeV/fm 3 )^ 2 , which gives us, from relation (|21|1 a 
value of ~ 10 18 Gauss. In this extreme case, the magnetic 
stress (B 2 /8n) ~ 10 35 dynes /cm 2 that is smaller but at 
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the same order of the matter pressure and can contribute 
to the stellar mass. 

In normal matter, a non-potential magnetic 
field is held inside a conductor in the form of 
stresses in the lattice. There the critical stress 
is smaller than unity, and consequently B 2 would 
be much less than the pressure P. However, in a 
highly dense and a complicated system like a neu- 
tron star where forces are close to the strong in- 
teraction regime, the limit of the neutron matter 
to support steady B field is still unknown and is 
expected to be many orders of magnitude higher 
than that of normal matter. 

The critical field as a function of the depth of the star 
is shown in Fig. @. The choice of different EOSs will 
change the number coefficients, but the overall order of 
magnitude will remain the same. It is also shown in Fig. 
Pjl that even for the maximum depth of the pressure gra- 
dient, the magnetic field is calculated to be still of that 
order of magnitude. So, the critical magnetic field depends 
essentially on the intrinsic scales of the neutron star pa- 
rameters, such as mass, pressure, radius, etc. 



1.8 - 

[ Emission of Bursts 

1.6 - I 
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Radius [km] 



Fig. 2. Critical magnetic field (B) line defined by the con- 
dition dp r /dr — and given by Eq. J5DJ. 



It is worth mentioning at this point, that the 
strong B field on dp/dr (of Ea.(|18|0 at r indicates 
that an entire shell of volume 47rr 2 dr should have a 
magnetic field stronger than the critical field limit. 
However, alignment of the magnetic moments of 
the nucleons, which is the possible reason for the 
formation of the strong field, is not guaranteed 
throughout the entire shell at the same time. So, 
the small pockets inside the shell which gets mag- 
netised, are referred here as the magnetic bubbles, 
and the dp/dr of Ea. (|18|) acts only on these small 
volume elements. The magnetic fields in these bub- 
bles are local. 



4. Gamma-Ray emission process 

One possible consequence of this general relativistic ef- 
fect coming from strong magnetic fields can be a mecha- 
nism for the emission of gamma-ray flares in SGRs. Energy 
emitted during each burst process in a SGR is of the order 
of 10 36 to 10 40 erg with some special exceptions of about 
10 44 erg (the March 5, 1979 event in SGR 0526-66 and 
the August 27, 1998 event in SGR 1900+14). 

Fig. J5J indicates the region below the curve as the 
stable region, where there is no reversal of sign in the 
modified TOV. Thus, this is the limit of the magnetic 
fields in the stars having no burst. When the field limits 
exceeds the critical value (as shown in Fig. J5J) of the 
magnetic field, then we have flares of gamma-rays from 
the stars. 

How much of the stellar matter is needed to produce a 
burst of the order of 10 36 erg? Typically, if we consider a 
small bubble of 1 cm radius, then the volume of the bubble 
is 10 -18 times than that of a neutron star of radius 10 km. 
For a bubble of radius 1 cm, the total energy released will 
be 

£ = r Z— ~ (1)3 x (io 18 ) 2 erff ~ 10 36 er<? (22) 
6 

where we have used the scale of the order of 10 18 Gauss 
as obtained for the critical field strength. 

So, to match with the energy released by the SGRs (~ 
10 36 erg/s), the bubble should have a volume of ~ 1 cm 3 . 
There is a little variation in the scale due to the position 
of the bubble inside the star. Another variation comes 
from the choice of the EOS of the matter inside the star. 
These collected variations however changes the numbers 
by maximally one order of magnitude, but the general 
features will be restored. We have considered a simple 
model where magnetic field condense in a small 
sphere inside the neutron star. So, this gives a very 
general picture in our model of freezing of magnetic field 
inside small droplets or bubbles, negligibly smaller than 
the volume of the star, and emission of gamma-ray flares in 
the SGRs, from the effect of relativistic equations. It can 
be equally applicable to any and every neutron/strange 
star model having a strong magnetic field. 

5. Conclusions 

We have discussed a general relativistic effect, due to 
strong magnetic fields, that can induce a radial instability 
inside magnetars, with the subsequent release of gamma 
rays as an outcome. We have estimated the strength of the 
magnetic field in order to drive such an effect. The high 
magnetic field that is necessary to change the chemical 
potential of the particles in the nuclear matter of the neu- 
tron star has been calculated to be more than 10 18 Gauss 
and the maximum mass of the entire star, due to this high 
field, increases by a very small percentage (e.g., see Table 
2 & 4 of Cardall et al., 2001). So, for our model, a field 
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of the order of 10 tolO 18 Gauss inside a small bubble of 
radius 1 cm, which is more likely to make a burst, the in- 
crease of chemical potential inside the bubble due to very 
strong held, will leave very little effect on the chemical 
potential in the rest of the matter of the star where the 
held is comparatively weak. It is also worth emphasising 
the fact that a freezing of the magnetic held in the small 
volumes, is more likely due to alignment of magnetic mo- 
ments of the nucleons, rather than due to any convective 
mechanism (e.g., Table 4 of Cardall et al., 2001), and can 
be continuously produced. They do not need to form just 
during the birth of the neutron star and survived for thou- 
sands of years. Besides, a burst of this kind in the interior 
of a compact star may produce a starquake as those al- 
ready observed in the SGRs. This relativistic effect due 
to the anisotropic nature of the stress-energy momentum 
tensor inside a strongly magnetised NSs or magnetars, as 
the progenitors for the central engine of the giant hares in 
SGRs seems to having been overlooked in previous studies. 
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